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for the high temperature superconductors has been done. How far the em- 



ployment of an asymptotic expression for the density of states underestimates 

s 

1 ■ the thermodynamic parameters, as for example, the critical temperature, the 

gap at zero temperature, and the known universal ratios of the BCS model, 



is discussed. Analytical expressions obtained for some thermodynamic func- 



^ . tions are compared with the numerical results for which the exact 2D density 

of states has been considered. Approximate analytical expressions for the 
temperature dependence of the reduced gap, the superconducting electronic 
specific heat, and the critical magnetic field have been obtained and compared 
with the corresponding numerical results. The validity of the different approx- 
imations and of the equations obtained for the thermodynamic functions in 
the framework of the van Hove model are also analyzed. 
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I. INTRODUCTION 



For electron-phonon (e-ph) superconductors, the Eliashberg equations (EE)0 allow the 
exact calculation of the temperature dependence of the thermodynamic functions. To solve 
EE one needs to know the e-ph spectral distribution function, a 2 F(u), and the coulomb 
repulsion pseudopotential parameter, fi*. The function a 2 F(uj) has been obtained from the 
inversion of tunneling data and calculated theoretically for several e-ph superconductors. 
As a rule, a very good agreement has been obtained between experimental and theoretical 
calculations. The parameter fi* is more difficult to deal with in a precise way. The common 
practice is to fit it to the experimental critical temperature, T c , through the linearized EE, 
valid at T c . These data are then used to solve the non-linearized EE from which the free- 
energy difference as a function of the temperature, T, follows!!. 

BCS theory is the first solution found that gave the clue to the explanation of the e- 
ph superconductivity. It is a weak coupling limit that, as it is to be expected, deviates 
sometimes strongly, from the experimental results even for some medium e-ph coupling 
superconductors^. An important role, as reference values, has been played by the universal 



ratios (UR), R x = 2A(0)/k B T c , R 2 = AC(T c )/C en (T c ), and R 3 = H c (0)/y/N(0)A(0), where 
A(T) is the gap function, AC(T C ) = C es (T c ) — C en (T c ) is the jump of the electronic specific 
heat at T c , given by the difference between its value in the superconducting state C es (T c ) 
and in the normal state C en (T c ), H c is the thermodynamic critical magnetic field, and N(0) 
is the density of states (DOS) at the Fermi energy. For e-ph conventional superconductors, 
the BCS UR do not depend on any parameter: R\ = 3.52, R 2 = 1.43, and R3 = l^/li. The 
BCS theory was mainly used as a reference, as a limit that helped characterize the relative 
strength of the e-ph coupling for conventional superconductors (CS). 

In the new high-Tc superconductors (HTS), the situation is quiet different since the 
mechanism driving the superconducting phase transition is still unknown. Actually, it is 
amazing how much it has been established about HTS without knowing the mechanism. 
Two main attitudes have been adopted^: to assume a mechanism and to calculate the ther- 
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modynamics that follows, or to incorporate some known facts into BCS theory which makes 
no assumption on the specific mechanism. Several known facts have been discussed within 
the BCS framework as, for example, the symmetry of the gap and the influence of the dimen- 
sionality on the phenomenon of superconductivity. In this sense, the van Hove BCS (v-BCS) 
scenario has been formulated by different authorsi'0. Both, analytical and numerical treat- 
ments have been carried out leading to results that present noticeable differences in the 
thermodynamic functions. Tunneling experiments and the dependence of the specific heat 
on temperature@, have been analyzed following the BCS formalism. For example, the con- 
tributions that come from the Cu02 planes, CuO chains, and c-axis, have been decoupled^. 
The one coming from the Cu02 plane seems to be the most important one. Most of the 
work has been devoted to it. In this sense, it is necessary to understand very well super- 
conductivity on the Cu02 plane and to analyze in detail how the different approximations 
used describe the phenomenon. In other words, it is important to be aware whether certain 
deviations from the experimental results come from physics or from mathematics. 

In this paper we want to analyze three types of solutions to the van Hove scenario, i.e., 
the analytic ones and their approximations, the numerical "asymptotic" solution that makes 
use of the asymptotic behavior of the two-dimensional density of states, and what we will 
call the "exact" solution where the full elliptic integral of the first kind, taking place in the 
2D DOS, is considered and treated numerically. The calculations presented here assume an 
s-wave symmetry for the gap in the Cu02 plane since it is the most widely studied case, 
although it is known not to be the experimental verified result. Nevertheless, this is not 
important for the sake of our purpose and the analysis obtained does not change in any 
essential way whether we introduced or not a d-wave or a mixing of s and d gap symmetries. 

The rest of the paper is organized as follows. In sec. 2, we compare the analytic, the 
asymptotic and the exact solution within the van Hove scenario for the gap, the critical 
temperature, the temperature dependence of the gap, the electronic specific heat, and the 
thermodynamic critical field. Discussions and conclusions are considered in sec. 3. 
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II. BCS MODEL WITH A VAN HOVE SINGULARITY IN THE DOS 



The layered structure of the new superconducting materials leads to almost dispersionless 
electronic spectrum on the c-axis@ and therefore, as a first approximation, it can be modeled 
as a stack of decoupled CuC>2 planes. This constitutes the so called van Hove scenario. For 
the CuC>2 planes, we take the following electronic dispersion relation: 

= — 2t(cos k x a + cos k y a) + At' cos(k x a) cos(k y a) , (1) 

where t and t' are the coupling parameters between the Cu-0 atoms and 0-0 atoms respec- 
tively, k x , k y are wave vector components, and a is the lattice constant. The origin of the 
energy is at the Fermi level, Ef- The resulting topology of the electronic band structure in 
the first Brillouin zone for t' = (a) and t' = O.At (b) are shown in Fig. 1. The saddle points 
at k = (tt /a, 0) and (0, n/a) are clearly seen yielding to the van Hove singularity (vHs) which 
enhances the DOS. The parameter t' ranges between and 0.5 t and it is responsible for 
the observed changes in the shape of the Fermi surface (compare Fig. 1(a) to Fig. 1(b)). 
It is worth noticing that for a fixed carrier concentration, the parameter t' drives the rela- 
tive position of the singularity with respect to Ef- Later on, we will show how almost all 
superconducting properties depend on this relative position. 
The single spin DOS for a stack of layers0 is equal to 



N{e) = ^_i= K 



16t 2 - (e - At 



A2 



(2) 



\\ 16t 2 (l + f) 

where K{m) is the elliptic integral of the first kind0. S is the area of the sample, n p is the 
number of layers, and iVo = \j2-n 2 t. The asymptotic form for Kim) in the neighborhood of 
the singularity with t' = is@ 

„ T/ s Sn v N , / 16 \ 

The equation above allows us to obtain analytical expressions for some of the thermodynamic 
functions. The BCS theory equation for the gap is 
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where Ep = + |A^,(T)| 2 , and Vg£, is the effective electron-electron interaction matrix 
element. In agreement with the BCS formulation we use the following parametrization for 



y i ; if N> \ £ %>\ < £ c 

; otherwise 



and e c is an energy cut-off. Following the standard BCS procedure, we get 



2 / tanh (Je 2 + A 2 (T)/2k B T 

- = de N(e) W y 

Vi J Je 2 + A 2 (T) 



(4) 



A. The gap and the critical temperature 

The numerical solutions of eq. (§) with a DOS given by eq. (HP will be referred to as the 
"exact" solution. When eq. (|) is used instead, we will talk about "asymptotic solution". 
We, first, compare the exact to the asymptotic DOS in the calculation of T c and A(0), in 
particular, their dependence on the dimensionless interaction parameter A = NqV, with 
V = ViUpS/a 2 . We present these results in Fig. 2. Here, a = 2t'/t, twice the 0-0 second 
nearest neighbors hopping parameter in units of the Cu-0 first nearest neighbors one. In 
Fig. 2(a), we present the dependence on A of the critical temperature T c and in Fig. 2(b), 
the one of the gap A(0), normalized to the cut-off energy e c . The curves correspond to 
a = 0.01, 0.04, and 0.07. Both parameters, T c and A(0), increase with A almost in the same 
way. For lower values of a (Cu-0 interaction stronger than the 0-0 one), T c and A(0) both 
are a very steep function of A. Hence, for a single plane, the dependence of the critical 
temperature and the gap on A is strong. It is worth noticing that T c of order 100K can 
be reached with A = 0.15 and a = 0.01. On the other hand, for higher values of a (a big 
relative hopping energy to the 0-0 atoms), i.e a = 0.07, A would have to increase as high 
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as 0.25 to obtain T c pa 100K. It is worth noticing that for a = 0.07, A needs to be of the 
order of 0.30 to get A(0)/e c « 1. 

In Fig. 3, we show the exact solution for the universal ratio, R\ = 2A(0)/ksT c as a 
function of A, for the same values of a considered above. At A=0.15 and a=0.01, i?i=3.82; 
while for higher values of A, R\ tends to 4 irrespective of the parameter a. Notice that, 
according to these results, the asymptotic value is reached for A's that might be very high 
compared to the expected ones in the cooper oxide sup er conduct or s0. On the other hand, 
at low A values, we get the usual BCS value, i.e., R\ = 3.52. On the overall, these results 
show that for any value of A and any chosen parameter a, 3.5 < R± < 4. Therefore, in any 
case, the 2D character of superconductivity in the cooper oxides rises R\ by at most 0.5 
within the BCS formulation. 

Approximate analytic solutions are obtained by inserting eq. (|3]) into eq. 



A(0H64 ( exp « +In .(^)- 1 | 



fc B T c =32 t exp(l-^ + ln*Gl)-\) ' (6) 
from which a parameter-independent universal ratio follows: 

The numerical calculations using the exact DOS converge to Rx = 4 for high values of A, 
so the analytic solution (7) , represents an upper bound for this ratio (see Fig. 3). The Ri 
value is smaller than those reported experimentally, where R\ ranging between 4 and 8 are 
foundlll. This disagreement is well outside the experimental error. 

In Fig. 4, the dependence of T c on the dimensionless interaction parameter, A, is studied. 
The solid lines are numerical solutions obtained from eq. (Q) with the DOS given by eq. (2) 
and t' = 0. The dot-dashed lines correspond to solutions using the asymptotic behavior of 
the DOS (eq. (3)). Three set of curves are presented in each case: e c = 20 meV (1), 35 meV 
(2), and 50 meV (3). As it can be seen, the cut-off energy has a strong influence on T c . At A = 
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0.20, for example, the solid lines range between WOK and 300i^ for 20meV < e c < 50meV. 
It is difficult to decide for a precise value of e c which actually acts as a free parameter and 
quantitative results are subject to a physical justification. In Fig. 4, we also show the results 
using the exact DOS (solid lines). Compare them to the asymptotic expression (dot-dashed 
lines). The exact DOS gives, for the same value of e c , always higher values for T c . At 
A=0.15, for example, the exactT c is about twice the one obtained from the asymptotic DOS 
given by eq. Furthermore, the analytic solution given by eq. (6) using e c = 50meV 

(dashed line) gives values for T c that are just (10%) above the numerical results using the 
asymptotic DOS but quite lower than the exact solution. So, the effect of the vHs is quite 
underestimated by using in either way the asymptotic form for K(m). Any analysis of HTS 
based on BCS theory should take into account these very strong dependences on several 
parameters before reaching any even qualitative conclusion. This is a main result of this 
paper. We will illustrate it further. 



B. Temperature dependence of A(T)/A(0) 

The temperature dependence of the reduced gap, A(T)/A(0), on the reduced tempera- 
ture T/T c is presented in Fig. 5. For the sake of comparison the Muhlschlegel modefl is 
also shown(in circles). It is found that the temperature dependence of A(T)/A(0) in the 
framework of the models considered here present the known universal character irrespective 
of the value of the parameters a, e c and A. 

In the following, approximate analytic expressions for the gap at low temperatures and 
close to T c are derived. In the low temperature regime and taking t' = 0, eq. (|]) can be cast 
in the form I\ — I2 — —2/3, with 

h= JdeN(e) 1 ; = / deN(e) I (j = 2,3) , 

where f{E) is the Fermi-Dirac distribution function. Using eq. (||D, the difference I\ — I2 is 
calculated in a similar way to that used in ref0: 
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h - I, « (1 + In 16) In + I f _ I f ^ 

V ; V A(0) / 2 \ At j 2 \ At j 



h oc 



(8) 



Collecting the results for Ii — I 2 and 1% and taking into account that R\ = A, the following 
analytical expression for the gap, A(T), is obtained 

A(T) . 3 / /ttT 10' 



1 



exp 



2T, 



K\ — dx In 



x 



x 



(9) 



A(0) * 4 \V T c ' 3 / ~' r V T 

Let us look now at the behavior of the gap near T c . For this purpose eq. (01) is rewritten as 
K l -K 2 = K 3 with 

16\ tanh(/3 c x/2) 



K 2 = I dx In ( — 

x J 



16\ tanh(#r/2) 



x 



K 3 = - W £ / dx In (- 

n=l X V X 



[(271 - 1) 2 7T 2 + (3 2 X 



212 



(10) 



where /3 C = At/k B T Cl x = e/At, y = A(T)/4t, u = e c /At, f3 = At/k B T, and the following 
identity has been used 

tanh - = Ab } rr— — . 

^2/ ^ (2ti - l) 2 vr 2 + 6 2 

Since if 3 decreases strongly for large values of /3x, the upper limit can be extended to infinity 
and 

2 



ln(2n-l) ln(16/3/vr) + 1 



(2tj - l) 3 (2n - l) 3 
The difference K\ — K 2 can be obtained using the same approach of ref.0 and results i 



(11) 



m 



Jfi-jr a =(l + lnl6) ^[^)-\ 



i„2 lb!E\ - w (kBT: 



2t 



2t 



(12) 



The above equations are the starting point to develop an explicit equation for the gap close 
to T c . Differentiating K\ — K 2 in eq. (p~2|) and K 3 from eq. (fTTJ) with respect to T results in 



dA 2 (T) 



dT 



57r 



'7((3) klTc 



-9.38k 2 B T c 



(13) 



which is the same as c-BCS result. To derive the above equation the following identities 
have been employedill 



E 



:i-2- k )((k) 



{ (2n - 1)* 



E 



ln(2n - i; 



2- k C(k) \n2-(l-2- h )C(k) 



(14) 



^ (2n - l) fc 

where is the Riemann Zeta function, and C'(^) its derivative. From eqs. (|TT| ) and (12) 
follow the functional dependence A 2 (T) = F(T). Irrespective to the form that the function 
F(T) has, for T near T c and due to the smallness of A 2 (T) , this can be expand in powers 
of (T — T c ), thus 



(T-T c 



From eq. ([T~5|) , A 2 (T) can be cast into 

A 2 (T) = 9.38 k 2 B T c 2 (l - | 
and taking into account eq. (7) it is obtained that 

A(T > - 1.53 4 fl~ 



(15) 



A(0) "™ V* T t 

The expressions for A(T) given by eqs. (|[) and (|15D, are plotted in Fig. 5 in dashed and dot- 
dashed lines, respectively. It can be seen that eq. (§) is valid in the range < T/T c < 0.4 
whereas eq. (|l~5|) in the range 0.95 < T/T c < 1. 



C. Electronic specific heat 

The fermion specific heat, C en , in the normal state is given by 
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C en (T) = yj-, J de N(e) e 2 f(e) [1 - /(e)] . (16) 



-4t-4t 



In Fig. 6(a) the C en (T)/T dependence on T is presented for different values of the parameter 
a and ep = 0. For a = 0, C en /T has a singular behavior as the temperature decreases. On 
the other hand, for a > 0.03 and very low temperatures, the electronic specific heat presents 
a linear behavior with T. These features are consistent with those reported in ref.H, where 
the usual constant coefficient in the linear term of the specific heat becomes a function of 
temperature. 

An approximate analytic expression for C en will be now derived. After inserting eq. @ 
in eq. flilf) we get 

'k B r 



C en (T) = ^-RN k B T 



1.73 - In 



At 



(17) 



where R is the universal gas constant, n p S/a is assumed to be equal to the Avogadro's 



number, and the specific heat is calculated per mole. The first term in eq. (|T7| ) is the 
usual contribution from constant DOS. The second one, TlnT, is due to the vHs. A similar 
dependence has been reported in. ref.0. 

The superconducting-state electronic specific heat can be written as 



CesiT) = ^r 2 j de N{e)f{E)[l - f(E)][E 2 - \t^±] . (18) 



2 f^^rr^ 1^ 2 (T)- 

For the derivative of the square of the gap as a function of T we use eq. (|13]). We have 
calculated the dependence of C es {T) /C en {T c ) on T/T c , for three different values of a, (e c = 20 
meV, A = 0.1). We get a universal behavior consistent with the one obtained for A(T)/A(0). 
Inserting the eq. @ into eq. ([18]) results in 

C es (T) = ARN k B T c ^2vr exp (-^) [s.46 + \ In (^) , 

where a value of k B T c /4t = 0.015 has been assumed. This approximated analytical equation 
is valid at low temperatures (0 < T/T c < 0.4) where the gap is almost independent of 
temperature (see Fig. 5). Finally, for the ratio C es (T)/C en (T c ) follows: 
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C es (T) 
C en (T c ) 



7.15a 



27^ 3 

TlT 



exp 



T 



1.4 + — In 
10 



T 



Notice that the expected exponential behavior has been obtained. 
The specific heat jump at T c is given by 

At-At' 



AC(T C 



1 dA 2 (T) 



k B T c dT 



deN(e)f(e)[l-f(e)} . 



(19) 



In Fig. 6(b) we show the dependence of R 2 = AC(T c )/C en (T c ) on a. For low values 
of a, i?2 approaches the asymptotic behavior of 1.95 as the cut-off energy increases. At 
intermediate values of a, Ri crosses the c-BCS value towards lower values. At a = |f we 
get a discontinuity (in disagreement with ref. 0). a = e c /2t corresponds to the vHs energy 
shift from the Fermi level (ef = 0). The jump observed in the figure is due to the vHs (see 
eq. (|3P). Finally the curve approaches asymptotically the c-BCS value for higher value of 
a. At higher values of a a crossover to the c-BCS value of 1.43 is reached. 

In order to estimate the maximal value that R2 can reach, we take a = 0. Then, from 
eqs. (§), (0), dT7|), and (0) we get 



AC(T C 



9.38 RN Q k B 



1.96 - 1 - In (Mk 
2 \ At 



and 



AC(T C 



2.85 



1.96- \\n{k B T c /At) 
1.73 - \n(k B T c /At) 



C e n(T c ) 

Assuming k B T c /At = 0.015, we get an upper limit for R 2 : 

AC(T C 



R-2 



C en ('/ (■; 



1.95 . 



This result is slightly higher than the conventional BCS model value. The difference repre- 
sents the influence of the vHs. Even thought the value of 1.95 follows the trend to higher 
values in the UR, this is not high enough to be in agreement with those reported from 
experimental estimations (2 < R2 < 4.6)!. 
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D. Critical Magnetic Field 



Finally, the thermodynamic critical magnetic field, H C (T), can be evaluated from the 
thermodynamic relationship 

= F n (T) - F S (T) , (20) 

where the F n (T)(F s (T)) is the free energy in the normal (superconducting) stateS We use 
eq. (2) in eq. (20) and get the reduced temperature dependence of the reduced thermody- 
namic critical magnetic field H C (T) / H c (0)thaX appears in Fig. 7 for two different values of 
e c and A at a = (see figure caption for details). Curve 4 is the parabolic behavior known 
for conventional superconductors. Notice the strong departure from the parabolic behavior. 

After some straightforward but slightly lengthy algebra, we arrive at the following ex- 
pression 



H 2 C (0) 



2tt 



4 , /A(0)\ 



- - In 

3 



V 64t ) 



N /n A 2 (0) 

where Q = a 2 d and d is the interplanar distance. From the above equation it follows that 
H c (0) has a similar behavior as A(0) as a function of e c , A and a. Further, using the 
inequality — In a; > I — x, for x < 1 and A(0)/64£ <C 1, a lower limit for the UR R% = H c (0) 



R 3 = , SdB > 2.16^ . 



/ J N /Sl A(0) is obtained 



'JVo/fiA(0) 

This result represents another UR, whose value is higher than the 3D BCS model 
(2-y/7r), fixing a minimum for this ratio . 

The temperature dependence of H C (T) near T c , can be obtained in the following way. 
Near T c , A(T) is a small and F S (T) can be expanded in powers of A 2 (T). To second order, 
we get: 

F S (T) - F n (T) = - ai (T c )A\T) , (21) 
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1 



n , k 



1 



2 ^ [(2n-l) 2 7T 2 +^ c 2 £?l 2 



(22) 



The sum over k in eq. 
eq. ([11]), we get 

H?(T) 
8tt 



can be transformed into an integral over the energy, and using 



0.29 



a 2 d 



2.84 - In 



k B T c s 
At 



From this expression it follows: 



HAT) ( T 



(23) 



Equation ( p3| ) has a slope of 0.83 which is smaller than the one in c-BCS (1.74). The 



approximate H C (T) expression for T near T c , of eq. (p3|) is valid in the range 0.75 < T/T c < 1. 



III. DISCUSSION AND CONCLUSIONS 

In the present paper the thermodynamics of the BCS model within the vHs scenario 
is carried out using both the exact DOS and its asymptotic behavior (analytically and 
numerically). We were interested in the intrinsic error that derives from the use of the 
asymptotic approximation to the van Hove DOS as compared to the exact one. For that 
purpose we assumed an s-wave symmetry for the gap in spite of the experimentally shown 
d or (d+s)-wave symmetry. For our purpose this does not make any difference. Our study 
helps ping point to specific differences with experiment that come from the mathematical 
approximations used rather than from physics. A first conclusion is that the widely spread 
use of the logarithmic asymptotic behavior near the singularity of the elliptic integral of 
the frst kind that is obtained for the DOS in the van Hove scenario is a very disputable 
approximation. It can lead to results that differ substantially from the ones obtained when 
the exact DOS is used. So most of the analytical expressions (obtained in this way) are to 
be used with great care, in general, even to draw from them qualitative conclusions. 

Furthermore, even using the exact DOS of the van Hove scenario, we still get (within 
v-BCS) a picture that does not always agree with experiment. 
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Numerical computations of the gap at zero temperature A(0) and T c as a function of the 
dimensionless interaction parameter, A, were carried out for different values of a = 2t'/t, 
with e c , the cut-off energy, fixed ( see Fig. 2). T c and A(0) increase with A, but as a is 
increased, with A fixed, a strong decrease in both T c and A(0) is observed. As a increases, 
the range of possible values of A that are consistent with the experimental results, increase 
towards those often found in CS. An analysis of Fig. 2 shows that for A > 0.15 the ratios 
kBT c /e c and A(0)/e c are higher than one. In CS, this fact is understood as a signature of 
lifetime effects!. Thus, for small values of a and A > 0.15, many body effects should play an 
important role . This is consistent with the explanation given for Csi to account for some 
high values found for R±. The universal ratio, R\ = 2A(0)/A;bT c , ranges between 3.5 and 4 
for any values of a and A, well below the experimental reported value. The highest values 
of A(0) and T c are obtained at a = 0. We have derived analytic expressions in this limit 
to examine further this point. The asymptotic approximation underestimate A(0) and T c 
as compared to the exact result. It is interesting to notice further that there is but a small 
difference using Eq. 3 whether the problem is solved analytically or numerically (only by 10 
% , see Fig. 4). 

The temperature dependence of the reduced gap A(T)/A(0) on the reduced temperature 
T/T c presents the known universal behavior irrespective of the value of a, e c and A. Analytic 
expressions can be obtained for the gap at < T < 0.4 (eq. (9)) as well as at 0.95 < T/T c < 1 
(eq. (15)). The presence of the vHs leads to higher values of A(0) and T c than those reported 
for the CS even for low values of the dimensionless interaction constant A. Furthermore, 
the magnitudes A(0), T c and R\ are strongly affected by the departure of the vHs from the 
£p. This is an obvious result since the number of states around Ef (that contribute to the 
superconductor state) decreases as the singularity moves far away from the Fermi level. 

We have further calculated the specific heat. Our results appear in Figs. 6(a) and 6(b). 
We get for C es (T)/T an exponential decay at low temperature (for a > 0) in contrast with 
that reported iniiS. It is possible to show, that the corresponding width in the transition 
region (T close to T c ) decreases as a increases in agreement with experimental evidences^. 
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On the other hand, in the range of temperatures close to T c , and due to thermally-created 
quasiparticle excitations, with roughly the same weight for all directions of k, qualitative 
differences with respect to c-BCS behavior are not found and the specific heat jump has a 
sharp maximum at T c . We get a linear dependence of AC on A(0) instead of the quadratic 
dependence reported In this region a broadened step-like discontinuity, as it has been 
found in some CS, is presented. For UR R2 = AC (T c ) / C en (T c ) we found that even though 
the value that we have obtained is higher than that of the c-BCS it is not enough to be in 
agreement with experimental reported that ranges between 2 and 4. 

The reduced thermodynamic critical magnetic field H C (T) / H c (0) as a function of T/T c 
shows a marked departure from the characteristic parabolic dependence observed in CS. 
Another important feature of the curve H C (T) / H c (0) is the change of curvature that takes 
place close to 0.3T C . The origin of this departure from the parabolic law is associated with 
the temperature dependence of the specific heat in the superconducting state. In ref.il a 
dependence like yT — (T/T c ) 2 is proposed instead of the parabolic law. Finally, a minimum 
value for the ratio R 3 = H c (0)/ ^Jn /QA(0) is obtained (3.83) which is higher than that 
reports in the c-BCS theory. 

Our two main conclusions are: First, the widely used logaritmic form of the DOS in 
the van Hove scenario is a very disputable approximation. Second, the van Hove scenario, 
even using the "exact" DOS although it enlightens quite a lot, might need more elements 
to account properly for the thermodynamics of HTS. It seems that other contributions 
neglected by the van Hove scenario play an important roleB 

Acknowledgments 

D. Q. would like to thank the financial support and the hospitality of the Department 
of Physics at CINVESTAV, Mexico, D.F. 

*Permanent address: Departamento de Fisica Teorica, Universidad de La Habana, 10400 La 
Habana, Cuba 



15 



REFERENCES 

^.M. Eliashberg, J. Exp. Theoret. Phys. 11, 966 (1960). 

2 J.P. Carbotte, Rev. Mod. Phys. 62, 1067 (1990). 

3 J.M. Daams, J.P. Carbotte, M. Ashraf, and R. Baquero, J. Low. Temp. Phys. 55, 1 (1984); 
R. Baquero and J.P. Carbotte, J. Low. Temp. Phys. 51, 135 (1983). 

4 J. Ruvalds, Supercond. Sci. Technol. 9, 905 (1996); A.S. Davidov, Phys. Rep. 190, 191 
(1990). 

5 J. Friedel, J. Phys. (Paris) 48, 1787 (1987); J. Phys. Condens. Matter 1, 7757 (1989). 

6 J. Labbe and J. Bok, Europhys. Lett. 3, 1225 (1987); J. Bok, Physica C 209, 107 (1993). 

7 C.C. Tsuei, D.M. Newns, C.C. Chi, and PC. Pattnaik, Phys. Rev. Lett. 65, 2724 (1990); 
D.M. Newns, C.C. Tsuei, PC. Pattnaik, and C.L. Kane, Comments Cond. Mat. Phys. 15, 
273 (1992). 

8 R.S. Markiewicz, J. Phys. Condens. Matter 2, 665 (1990); Physica C 168, 195 (1990). 
9 J.M. Getino, M. de Llano, and H. Rubio, Phys. Rev. B 48, 597 (1993). 
°A.G. Goicochea, Phys. Rev. B 49, 6864 (1994). 

1 R. Baquero, D. Quesada, and C. Trallero-Giner, Physica C 271, 122 (1996). 
2 R. Grassme and P. Seidel, J. Supercond. 9, 619 (1996). 

3 D.M. Newns, C.C. Tsuei, and PC. Pattnaik, Phys. Rev. B 52, 13611 (1995). 

4 M. Liu, D.Y. Xing, and Z.D. Wang, Phys. Rev. B 55, 3181 (1997). 

5 A.M. Cucolo, C. Noce, and A. Romano, Phys. Rev. B 53, 6764 (1996). 

6 W.E. Pickett, Rev. Mod. Phys. 61, 433 (1989); Z.X. Shen and D.S. Dessau, Phys. Rep. 
253, 1 (1995). 

16 



E.N. Economou, Green functions in quantum physics, Springer Verlag, Berlin (1983). 

M. Abramowits and I. Stegun, Handbook of Mathematical Functions, Dover Publications 
Inc. N.Y. (1970). 

T. Hirata and Y. Asada, J. Supercond. 4, 2 (1991). 
B. Muhlschlegel, Z. Phys. 155, 313 (1959). 

I.S. Gradsteyn and I.M. Ryzhik, Table of Integrals, Series and Products, Academic Press, 
San Diego (1980). 

A. Junod in: Physical properties of high temperature superconductors II, ed. by D.M. 
Ginsberg, p. 13, World Scientific, Singapore (1990). 

J.R. Schrieffer, Theory of superconductivity, W.A.Benjamin Inc. (1964). 

B. T. Geilikman, V.Z. Kresin, and N.F. Masharov, J. Low. Temp. Phys. 18, 241 (1975). 

N. Momono and M. Ido, Physica C 264, 311 (1996); N. Momono, M. Ido, T. Nakano, and 
M. Ota, Physica C 235-240, 1739 (1994). 

T.E. Mason, G. Aeppli, S.M. Hayden, A. P. Ramirez, and H.A. Mook, Phys. Rev. Lett. 
71, 919 (1993). 

S. Dorbolo, M. Houssa, and M. Ausloos, Physica C 267, 24 (1996). 
J.H. Cho, Z. Hao, and D.C. Johnston, Phys. Rev. B 46, 8679 (1992). 



17 



FIGURES 

FIG. 1. Topology electronic band structure in the first Brillouin zone following eq. (1) for t' = 
(a), and t' = OAt (b). The energy is given in eV and the calculations correspond to t = 0.25 eV. 

FIG. 2. Dependence of the critical temperature, T c (a) and the reduced gap at zero temperature, 
A(0) (b) on the dimensionless interaction constant A = NqV, at e c = 20 meV and a = 0.01 (solid 
line), 0.04 (dashed line) and 0.07 (dot-dashed line). 

FIG. 3. Dependence of the ratio R± = 2A(0)/fcsT c on the dimensionless interaction constant 
A = NqV, for the same values of a and e c of fig. 2. The c-BCS and v-BCS universal ratios are 
presented by thin lines. 

FIG. 4. Critical temperature T c as a function on the dimensionless interaction constant 
A = NqV, for e c : 20 meV (curves 1), 35 meV (curves 2), and 50 meV (curves 3). The calcu- 
lations with the exact DOS, eq. (|2|) and its asymptotic behavior, eq. (|3|) are plotted by solid 
and dot-dashed lines, respectively. The approximate expression given by eq. (6) is represented in 
dashed line for e c = 50 meV. The value of a = has been chosen for all calculated curves. 

FIG. 5. Dependence of the ratio A(T)/A(0) on the reduced temperature T/T c . The T — > 
(curve 1 by dashed lines) and T — > T c (curve 2 by dot-dashed lines) approaches according to eqs. 
@ and (|l^), respectively are presented (for details see text). For comparison the c-BCS resul1@ 
is depicted in solid circles. The calculation correspond to a = 0. 

FIG. 6. Electronic specific heat in the normal state C en in units of T as a function of T for 
different values of a (a); The jump ratio AC (a) / C en (T c ) as a function of a for different values of 
e c at A = 0.1 (b). For comparison the c-BCS and the v-BCS universal ratios are presented by thin 
lines. 
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FIG. 7. Temperature dependence of H c (T)/H c (0) on T/T c obtained from eq. (^) with the 
DOS given by eq. @ for e c =35 meV, A =0.1 (curve 1); e c =20 meV, A =0.1 (curve 2); e c =35 
meV, A =0.02 (curve 3). For comparison the empirical parabolic law (curve 4) of the c-BCS (see 
ref£j) is shown. 
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